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We consider the spin-averaged nucleon forward Compton scattering amplitude in heavy 
baryon chiral perturbation theory including all terms to order 0(q 4 ). The chiral predic- 



tion for the spin-averaged forward Compton scattering amplitude is in good agreement 
with the data for photon energies ui < 110 MeV. We also evaluate the nucleon electric 
and magnetic Compton polarizabilities to this order and discuss the uncertainties of the 



various counter terms entering the chiral expansion of these quantities. 
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I. INTRODUCTION AND SUMMARY 



Compton scattering off the nucleon at low energies offers important information 
about the structure of these particles in the non-perturbative regime of QCD. The spin- 
averaged forward scattering amplitude for real photons in the nucleon rest frame can be 
expanded as a power series in the photon energy uj, 

r / 2\ 2 4 ^ ' 

/i(o; ) = a + ai u + a 2 uj +... 

where e^j are the polarization vectors of the initial and final photon, respectively, and 
due to crossing symmetry only even powers of u occur. The Taylor coefficients encode 
the information about the nucleon structure. The first term in eq.(l), ao = — e 2 Z 2 /4irm, 
dominates as the photon energy approaches zero, it is only sensitive to the charge Z 
and the mass m of the particle the photon scatters off (the Thomson limit). The term 
quadratic in the energy is equal to the sum of the so-called electric (a) and magnetic 
(/?) Compton polarizabilities, a± = a + (3. Corrections of higher order in uj start out 
with the term proportional to a 2 - 

In ref.[l], we calculated the electric and magnetic polarizabilities of the proton 
and the neutron in heavy baryon chiral perturbation theory (CHPT) to order q 4 (with 
q denoting a small external momentum or meson mass). The chiral prediction for the 
electric polarizabilities a p and a n are in good agreement with the available experimental 
data. Furthermore, we found a non-analytic term of the type In M n (with M n the pion 
mass) which cancels parts of the large positive A(1232) contribution to the magnetic 
polarizability f3 p (in case of the neutron, the coefficient of the lnM^ term is somewhat 
too small). This is a novel effect which allows to understand the relative smallness of 
the nucleons' magnetic polarizabilities. The aim of the present paper is twofold. First, 
we wish to consider the spin-averaged forward Compton amplitude as a function of the 
photon energy to investigate the importance of the a; 4 , a; 6 , ... terms and to confront it 
with the available experimental data. Second, as already stressed in ref. [1], we have to 
present a detailed analysis of the uncertainties entering the theoretical predictions for 
the polarizabilities at order q 4 . 

Heavy baryon CHPT is a systematic expansion of low-energy QCD Green functions 
in small external momenta and quark masses. It was previously used in the calculation 
of the baryon mass spectrum [2], axial currents in flavor SU(3) [3], nuclear forces [4] and 
many other observables (for a review, see ref. [5]). A systematic analysis of QCD Green 
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functions for nucleons in flavor SU(2) can be found in [6]. We will heavily borrow from 
that paper. We will work in the one loop approximation [6] which still allows to include 
all terms of chiral order q 4 . The two loop contributions start at 0(q 5 ). The nucleon 
electromagnetic polarizabilities belong to the rare class of observables which to leading 
order in the chiral expansion are given as pure loop effects [7] but at order q 4 counter 
terms enter. Their coefficients have already been determined in [1] by either direct fits to 
existing data for other reactions or by making use of the resonance saturation hypothesis 
[8,9]. This will be discussed in more detail below. 

The pertinent results of this investigation can be summarized as follows: 

(i) The spin-averaged forward Compton amplitude for the proton is in agreement with 
the data up to photon energies u ~ M n . It is dominated by the quadratic contribu- 
tion, i.e. to order q 4 in the chiral expansion the terms of order u 4 (and higher) are 
small. Similar trends are found for the neutron with the exception of a too strong 
curvature at the origin. 

(ii) The electromagnetic polarizabilities are in good agreement with the data (with the 
exception of j3 n ), see eqs.(28) and (29). We have discussed the theoretical uncer- 
tainties to order q 4 and found that the electric and magnetic polarizabilities can be 
predicted with an accuracy of ±2 • 10 -4 fm 3 and ±4 • 10 -4 fm 3 , respectively. An 
improved theoretical prediction is only possible if one can pin down the coefficients 
of some counter terms more accurately. 

The paper is organized as follows. In section 2, we briefly discuss the effective 
meson-baryon lagrangian underlying our calculation. In section 3 we introduce the spin- 
averaged forward Compton tensor which contains the information about the amplitude 
fi(u> 2 ) as well as the electromagnetic polarizabilities a and /3. Section 3 contains the 
results for the spin-averaged forward Compton scattering amplitude. The theoretical 
predictions for the electromagnetic polarizabilities are then given in section 4. 

2. EFFECTIVE LAGRANGIAN 

The basic tool to investigate the low-energy behaviour of QCD Green functions is 
an effective Lagrangian formulated in terms of the asymptotically observed fields, the 
Goldstone bosons of spontaneous chiral symmetry breaking (pions) and the matter fields 
(the baryons). The construction principles for building a string of terms with increasing 
number of derivatives are chiral symmetry, lorentz invariance and the pertinent discrete 
symmetries of the strong interactions as well as the systematic chiral power counting 
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scheme (see [5] and references therein). For our purpose, it suffices to stress that we 
work in the one loop approximation including all terms of order q 4 . As stressed in [3], 
n loop graphs start to contribute at order q 2n+1 if one inserts vertices from the lowest 
order effective Lagrangian C^. Notice that we consider the two flavor case and work 
in the isospin limit m u = ma- The pion fields are collected in the SU(2) matrix 

U(x) = exp[zf ■ n(x)/F] = u 2 (x) (2) 

with F the pion decay constant in the chiral limit. The nucleons are considered as very 
heavy [2,3,4,5,6]. This allows for a projection onto velocity eigenstates and one therefore 
can eliminate the troublesome nucleon mass term (of the Dirac lagrangian for nucleons) 
thereby generating a string of vertices of increasing power in 1/m, with m the nucleon 
mass (in the chiral limit). 1 -* The effective Lagrangian to order 0(q 4 ), where q denotes a 
genuine small momentum or a meson (quark) mass, reads (we only exhibit those terms 
which are actually needed in the calculation of forward Compton amplitudes) [1] 

£[ 1} N =H(ivD + g A S-u)H 

+ c aU - U - i ^{S.D,v- U }-±[S*,S»]((l + K v )f+ 
+ ^(Sa n + 6p n - 5a p - 5P p )Hf+J^Hv^v x 

(3) 



2 

- ^(5cx n + 6P n )HHTr(f+j¥)v»v x 



with 



u u = iu^V JJv) 

(4) 

/+ = (dpA,, - d v A^)(uZu^ + u^Zu) 

where H denotes the heavy nucleon field of charge Z = (1 + ts)/2 and anomalous 
magnetic moment k = (k 8 + tsk v )/2, the four-velocity of H, the covariant spin- 
operator subject to the constraint v ■ S = 0, the covariant derivative acting on 



^ Later on, we will identify m with the physical nucleon mass. This is a consistent 
procedure to the order we are working. 
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the pions, _D M = + T M the chiral covariant derivative for nucleons and we adhere to 
the notations of ref.[6]. The superscripts (1,2,4) denote the chiral power. The lowest 
order effective Lagrangian is of order 0(q). The one loop contribution is suppressed 
with respect to the tree level by q 2 thus contributing at 0(q 3 ). In addition, there are 
one loop diagrams with exactly one insertion from These are of order q 4 . Finally, 

there are contact terms of order q 2 and q 4 with coefficients not fixed by chiral symmetry. 
For the case at hand C^ N does not have to be considered explicitely. To order q 4 its 
vertices enter only tree diagrams and it therefore contributes to the polynomial part 
of the amplitude. The most general polynomial piece at order q 4 is, however, already 
given by the contact vertices of C^ N . Notice that some coefficients in C^ N related to 
the *yyNN and ^ixNN vertices are fixed from the relativistic theory by the low-energy 
theorems for Compton scattering (ao = — e 2 Z 2 / Aum) and neutral pion photoproduction, 
respectively. This is discussed in some detail in ref.[6]. The unknown coefficients we 
have to determine are ci, C2 and C3 characterizing a higher derivative ttttNN vertex as 

well as the four low-energy constants 5a p , oa n , 5/3 p and d(3 n from C^ N . We have not 

(2) 

exhibited the standard meson Lagrangian C^i, 

To calculate all terms up-to-and-including order o 4 , we have to evaluate all one 
loop graphs with insertions from and those with exactly one insertion from C^ N . 
While the former scale as g 3 , the latter constitute the new contributions of order q 4 . 

Furthermore, there are the tree diagrams related to C^ N which are also new. It is 

(2) 

worth to stress that in the one loop diagrams involving the anomalous magnetic 
moment of the nucleon appears, since the photon-nucleon vertex is expanded in the 
external momentum. We will come back to this point later on. Finally, it is mandatory 
to expand the leading order O(q) effective vertices as well as the nucleon propagator 
to include all relativistic corrections of order 1/m. The resulting Feynman rules for 
the pertinent vertices and propagators are shown in fig.l. Note that we have left out 
all vertices which give zero contribution either due to spin-averaging or due to isospin 
algebra. 

In eq.(3) we have expressed the effective lagrangian in terms of the physical param- 
eters m, k s , k v and so on. For the loop diagrams this is legitimate to the order we are 
working. Concerning tree diagrams which are responsible for the Thomson amplitude 
ao = — e 2 Z 2 / Aiim, the term proportional to c\ will shift the nucleon mass in the chiral 
limit by — Ac\M 2 to the physical nucleon mass at this order. 
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3. SPIN-AVERAGED FORWARD COMPTON TENSOR 



The object one has to study in order to determine the Compton amplitude fi(uj 2 ) 
and the electromagnetic polarizabilities is the spin-averaged Compton tensor in forward 
direction G, 



,2 



e^ = ^Tr[(l + y)T Ml >,A;)] 

= A(u) g^ v + B{u) k^k v + C(u) {k^v v + v^k v ) + D{u) v^v v 

where v and k denote the nucleon four-velocity (v 2 = 1) and the photon momentum 
(k 2 = 0) and to = v ■ k. T^{v, k) is the Fourier-transformed nucleon matrix element of 
two time-ordered electromagnetic currents, 

T^(v, k) = J d 4 xe ik - x < N(v)\TJ e ™(x) J?»(0)\N(v) > (6) 

Gauge invariance Q^ v k v = implies C(u) = —A{uj)/uj and D(u) = and therefore 
all information is contained in the two crosssing even functions A(co) = A(— uj) and 
B(u) = B(—u). This fact allows us to choose the "Coulomb gauge" e • v = for the 
polarization vector e of the photon when calculating the auxiliary quantity t v = 

A(u) e 2 + B(u>) (e • k) 2 . The gauge e • v = is very convenient since it leads to a drastic 
simplification in our loop calculation. In this gauge the leading order photon nucleon 
vertex vanishes and therefore many diagrams. The function A(u) is directly proportional 
to the spin-averaged forward Compton amplitude fi(uj 2 ) introduced in eq.(l), namely, 

A{u) = -4tt h{u 2 ) (7) 

The second function B(u) is not a physical amplitude, it only serves to calculate the 
magnetic polarizability. 

The electric and magnetic Compton polarizabilities are defined as follows: 

■= A"(0) 5 B(0) 

where the prime denotes differentiation with respect to uj. 

Our task is to calulate all contributions to A(u) and B(u) coming from CHPT 
at order q 4 . These are all one loop contributions with vertices from and those 

(2) 

with exactly one vertex fron C^^ as well as the most general polynomial contribu- 
tion at this order. The latter has the form A(u) po1 = e 2 Z 2 /m — 4ir(5a + 8(5) uj 2 and 
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B(u)) po1 = —47r5/3 where the constant term in A{ui) po1 is fixed by the low energy theorem 
for Compton scattering. Later we will determine the polynomial coefficients Sot, 8/3 by 
identifying them with the resonance contributions to the nucleon electromagnetic po- 
larizabilities. In making this identification we have excluded any contributions from the 
nucleon pole graphs to the polarizabilities. This prescription is in accordance with the 
analysis of the Compton scattering experiments where terms coming from the anoma- 
lous magnetic moment (nucleon pole graph) are explicitely separated from the ones of 
the polarizabilities by the use of the Powell cross section formula [15]. 

4. FORWARD COMPTON AMPLITUDE 

We wish to discuss the spin-averaged forward Compton amplitude A{uj) = 
—4irfi(u) for real photons (k 2 = 0) scattering off protons or neutrons. Making use 
of a once-subtracted dispersion relation and the optical theorem allows to express A{uj) 
in terms of the total photo-nucleon absorption cross section a tQt (uj) as 

e 2 Z 2 2u 2 f 00 j ,°tot(«/) 



n A , , e 2 Z 2 2u 2 „ f°° , ,<7 t 
ReA(u) = V duo — 



2 - u 2 (9) 
lmA(u) = -u>a tot (u) 

for the proton and the neutron. The threshold energy for single pion photoproduction 
is given by 

u = M 7r (l + M 7r /2m) . (10) 

Working in the Coulomb gauge e • v = 0, one has to calculate 9 irreducible one 
loop graphs with insertions from and 67 one loop diagrams with one insertion from 
C^jj. Furthermore there are the contact term graphs stemming from C^ N , which give 
rise to the most general polynomial contribution at order to A(oj) and B(u), namely, 
A(u) po1 = e 2 Z 2 /m - 47i(5a + SP)u 2 , B(u) po1 = -AnSp. 

Due to the simple structure of the loop integrals in the heavy mass formulation, 
one can give A P)n {uj) in closed analytic form. For the proton it reads 

e 2 

AJuj) = 47r(a + (3) p co 2 

m 

e 2 g 2 A M 7T f 3 1 . 1. r 1 11 2 ] 

+ i * { -- - -2 + (1 + -2 ) V 1 - z 2 + - arcsin^ + — z 2 } 
8ttF 2 { 2 z 2 v z 2 ' z 24 J 



e 2 g 2 A M 2 / 5 



z 



2 



+ 8^fi-2 + T (3K ' + 11) + 



2 



- (6 + k s )z + 



z 1 1 — z 



\/l — z 2 arcsin z 



1 ( 1 \ . 2 1 

+ - - Ks) arcsm z > 



(11) 
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with z = u/M^. Similarly, one has for the neutron 

e 2 g 2 A M n 



A n (u) =-Aix{a + f3) n u 2 + 



87TF2 



+ 



e 2 g\Ml\l z 2 (n n . 



( 3 1 , 1 N / ~ 1 

(« 



11 , 

arcsin z H -z" 

24 



2)^ + 



^ arcsin z 



+ - 72") 



arcsin z 



(12) 

Notice that the expressions for A Pjn (oj) diverge at oj = M x . This is an artefact of the 
heavy mass expansion. The realistic branch point coincides with the opening of the 
one-pion channel as given above. To cure this, let us introduce the variable 



c 



l + M n /2m M 7r (l + M 7T /2m) 



UJ 



(13) 



If one now rewrites A Pjn (uj) in terms of £, the branch point sits at its proper location 
and A Pin (C = 1) is finite. We have 



e 2 g 2 A M 7T f 3 1 , 1 N , 7 1 _ 11 - 



+ 



87r 2 mF2 \ 



+ 7T 



1 



1 



10 ^2 



11 



C ; 1 - C 2 



a/1 — C 2 arcsin £ 



— r arcsin C H C 

C 2 2C 24 s 



(i-C 2 ) 2 + i 

2C 2 V / W 2 ~ 
1 , 1 



± 



+ (i - (2 + K ,)C)v / l-C 2 arc S inC + i( f2 



.) arcsin 2 £ 



13, 



(14) 



where the '+/-' sign refers to the proton/neutron, respectively. The proper analytic 
continuation above the branch point ( = 1 is obtained through the substitutions 
a/1 — C 2 = — ^a/C 2 — 1 and arcsin £ = n/2 + zln(£ + a/C 2 — 1). The expressions given 
in eq.(14) differ from the ones in eqs.(ll,12) only by terms of order q 5 (and higher) and 
are thus equivalent to the order we are working. We should stress that in the relativistic 
formulation of baryon CHPT such problems concerning the branch point do not occur 
[10,11]. In the heavy mass formulation we encountered this problem since the branch 
point ujq itself has an expansion in 1/m and is thus different in CHPT at order g 3 and 
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We now present our numerical results. We always use the Goldberger-Treiman 
relation to express qaIF^ as g^n /m , with g n N = 13.40 the strong pion-nucleon coupling 
constant. We also use k s = -0.12, m = 938.27 MeV and e 2 /4n = 1/137.036. In 
fig. 2a, we show the real part of A p (u) normalized to the Thomson limit, e 2 /m = 3.02 
/ubGeV -1 , for the z and £ expansions in comparison to the data [12] for the central 
value of (a + 0) p = 14.0 • 10" 4 fm 3 of ref.[l]. Up to u ~ 100 MeV, the agreement 
of the prediction with the data is good. The corrections of order u 4 (and higher) are 
fairly small as shown in fig. 2b. One also recognizes the unphysical behaviour of the 
^-expansion around u> = M n . Similar statements hold for the neutron exhibited in 
fig. 3 using the value (a + f3) n = 21.2 • 10~ 4 fm 3 of ref.[l]. Using the deuteron data 
of Armstrong et al.[13], one finds at ui = 100 MeV, A n /A p = —0.35 compared to the 
theoretical prediction of —0.47. The difference is mostly due to the too large sum of the 
electric and magnetic polarizabilities of the neutron. 

For the proton, we have also calculated the real and imaginary parts of A(u) for 
C > 1. The imaginary part starts out negative as it should but becomes positive at 
ui ~ 180 MeV. This is due to the truncation of the chiral expansion and can only be 
overcome by a more accurate higher order calculation. Consequently, the real part 
(normalized to the Thomson limit) stays rather flat after the branch point as shown in 
fig.4. 

5. ELECTROMAGNETIC POLARIZABILITIES 

In [1] , we derived the following formulae for the electric and magnetic polarizabilities 
of the proton and the neutron (i = p, n) 



$Cg A 
4M W 



A C 
T 7T 



- c 2 ) In ^ + \ - 6c 2 + c+) 

v m 1 X 4 V 2m " ' _ 

(2£ki_ C2)1I1 ^ + i ( M9i + 2C2 _ <;+) 

v m ' A 4 v m ~ ' 



+ 5«[(A), 



(15) 



with 



C = 



= 4.36 • 10" 4 fm 2 , 



96tt 2 F2 

x n = 3 , y p = 71 



Vr 



X 



39 
37 

y 



p — 3 + k s 7 x n — 1 k s , y p — — + 6k s , y n — — 6k s , 

-801+402 + 403-1^- • 
2m 



(16) 
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Here, A is the scale introduced in dimensional regularization. The physical a.i and 
are of course independent of this scale since the renormalized counter terms Sa^(X) and 
8Pl(\) cancel the logarithmic A-dependence of the loop contribution. The corresponding 
renormalization prescription reads: 

*«* = ^2 U ~ Xj Pj + K(A), m = ^r 2 U- ^) + 6R(\) (17) 



7T 



m 



with 

_ X d ~ 4 
~ 16tt 2 



1 

-(7s — 1 — ln47r) 



(18) 



d-4 2 

The first term on the r.h.s. of eq.(15) is, of course, the result at order g 3 [6,7,10]. The 
results shown in eqs.(15) have the following structure. Besides the leading 1/M n term 
[6,7,10], 0(q 4 ) contributions from the loops have a \nM n and a constant piece ~ M°. As 
a check one can recover the coefficient of the In M n term form the relativistic calculation 
[1,2] if one sets the new low energy constants c« and k 8jV = 0. In that case only the 
1/m corrections of the relativistic Dirac formulation are treated and one necessarily 
reproduces the corresponding non-analytic (logarithmic) term of this approach. The 
term proportional to C2 lnM^ in eqs.(15) represents the effect of (pion) loops with 
intermediate A(1232) states [14] consistently truncated at order q 4 . We should stress 
that the decomposition of the loop and counter term pieces at order q 4 has, of course, 
no deeper physical meaning but will serve us to separate the uncertainties stemming 
from the coefficients accompanying the various contact terms. Notice that from now on 
we will omit the superscript V on 5a^(X) and 501 (X) appearing in eqs.(15). 

The numerical results for the various polarizabilities now depend on the knowledge 
of the contact terms C2, c + ,dai(X) and 5f3i(X) for a given choice of the scale A. As 
already stated, ideally the A-dependence from the loops is cancelled by the one from the 
corresponding contact terms, as detailed in [16], when one is able to fit all low-energy 
constants from phenomenology. We are not in that fortunate position but in some 
cases must resort to the resonance saturation hypothesis [8,9] to estimate some of the 
constants. In that case, the actual value of a certain low-energy constant is given as a 
sum of a dominating resonance contribution at the scale A equal to some resonance scale. 
This problem will eventually be cured when sufficiently many accurate low-energy data 
in the baryon sector will be available. Clearly, the contact terms we are dealing with 
fall into two categories. While C2 and c + only enter via higher order (in 1/m) vertices 
in the loop diagrams, ddi(X) and SPi(X) are "genuine" new contact terms of order q 4 . 
We will therefore discuss the sensitivity on these separately. Note, however, that only 
the total result of order q 4 is of physical relevance. 
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The constant c + can be related to the isospin-even S-wave pion-nucleon scattering 
length a+ [17], 

32 7 rM 2 F 2 (l-M 7r /m) V ; 

for a+ = -0.83 ± 0.38 • lO^/M^ with = 139.57 MeV the (charged) pion mass. 
Inserting the value of c + as given from eqs.(19) into eqs. (15,16), one finds that the 
contributions proportional to c + to the electromagnetic polarizabilities are negligible 
since they are less then 0.15 • 10" 4 fm 3 in magnitude. If we take only the leading order 
relation /87rM 2 , which is allowed to the order we are working, then the c + 

contribution to the electromagnetic polarizabilities is even smaller and only 0.06 • 10 -4 
fm 3 in magnitude. To determine C2, we make use of the resonance saturation hypothesis 
and get contributions from the A(1232) as well as the Roper AT* (1440), 



c A = 9A m 

° 2 4m 2 A 



n* Rg A m 

Co — 



m A + m 2 
ttia — m 

2 , 



(20) 



8(771^, — m 2 ) 



where Z parametrizes the off-shell behaviour of the spin-3/2 Rarita-Schwinger field in 
the AiW-vertex and we have used g^NA = ^g-n-N/V^ as well as g^NN* = ^fRg-irN /2. 
Empirically, Z is not known very accurately, —0.8 < Z < 0.3 [18], and R varies from 
0.25 to l. 2 ) Notice that the Roper contribution is more than one order of magnitude 
smaller than the one from the A. If one treats the A non-relativistically (isobar model), 
one arrives at 

for gA — 1.328 by making use of the Goldberger-Treiman relation. Altogether, c<i varies 
between 0.4 and 0.6 fm. The resulting polarizabilities depend only weakly on the actual 
value of C2 as shown in fig. 5. (for A = 1.232 GeV). As discussed before, there is some 
spurious sensitivity on the value of A as shown in fig. 6. The corresponding bands refer 
to the choice of g A = 1-26 or g A = 1.328 which are equivalent to the order we are 
working. While the neutron polarizabilities are quite insensitive to the choice of A, the 



2 - ) In ref.[19], a narrower range for Z is given under the assumption of g<i = which 
seems to be excluded by the data as stressed in [18]. We therefore use the wider range of 
Z given in [18]. 
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much larger coefficients x p and x' p induce some scale dependence for the proton case. 
Since most of the counterterms are in fact given by A exchange, we have chosen in [1] 
A = niA (which gives our best values). 

The A(1232) enters prominently in the determination of the four low-energy con- 
stants from C^Ij. Therefore, we will determine these coefficients at the scale A = m A . 
In particular, one gets a sizeable contribution to the magnetic polarizabilties due to 
the strong NA Ml transition. A crude estimate of this has been given in ref.[20] by 
integrating the Ml part of the total photoproduction cross section for single pion pho- 
toproduction over the resonance region, 

S^(m A ) = ^ [ ^a M1 M=7.0-10- 4 fm 3 (22) 

However, this number is afflicted with a large uncertainty. If one simply uses the Born 
diagrams with an intermediate point-like A, one finds 

= T^{ m KZZ: m2 - ^(2Y + 1) + MY + D] } (23) 

with gi the strength of the 7./VA coupling and the off-shell parameter Y is related to the 
electromagnetic interaction C^ NA (see ref.[18] for more details on this). These parame- 
ters are not very well determined, a best fit to multipole data for pion photoproduction 
leads to 3.94 < g x < 5.30 and -0.75 < Y < 1.67 [18]. The recent PDG tables give 
3.5 < <7i < 7.5 [21]. For comparison, SU(4) and chiral soliton models give g\ = k* = 5.0 
which was used e.g. in [10] together with Y = —1/4 to estimate the A contribution at 
order q 4 . Inspection of eq.(23) reveals that the large positive values of Y lead to very 
large negative contributions from the A in plain contradiction to the dispersive estimate 
of eq.(22). If one, however, assumes an universal off-shell parameter for the strong and 
electromagnetic iVA transitions, it is plausible to set —0.8 < Y < 0.3. In that case, 
5(3 A is always positive and varies between 14. and 0.5 as shown in fig. 7. A conservative 
estimate therefore is 

SPp(m A ) = SP^(m A ) ~ (7.0 ± 7.0) • 10" 4 fm 3 (24) 

invoking isospin symmetry. Clearly, the large range in the value for 8/3 A is unsatisfactory 
and induces a major uncertainty in the determination of the corresponding countert- 
erms. In [1], we choose the central value of eq.(24) as our best determination. In 
the framework of a non-relativistic calculation (isobar model), the corresponding Born 
term which generates 5(3 A does not depend on any off-shell parameter and one finds 
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e 2 # 2 /(187rm 2 (m A - m)) = 12 • 10" 4 fm 3 (for g x = 5). In ref.[22], the A(1232) was 
included in the effective field theory as a dynamical degree of freedom and treated non- 
relativistically (like the nucleon). It has been argued in [22] that the A Born graphs 
have to be calculated at the off-shell point ui = 0. This effect can reduce the large 
5f3 A by almost an order of magnitude. This is reminiscent of the off-shell dependence 
discussed before. It was already pointed out in ref.[10] that a relativistic treatment of 
the A(1232) also induces a finite electric polarizability at order q 4 . This contribution 
depends strongly on the 7A./V couplings gi and g 2 as well as the two off-shell parameters 
IF, 



8a A (m A ) 



18nm 2 m 



9\ 



+ 9W2 



92 

4 



mA + m 
m A {m — mA) 



+ 4r(m A (l + 2Y) -mY)) 



2(mA + Tn 
4m A + m A m + m 2 



+ (X + Y + AXY)m A - Y(l + 2X)m 
+ X(1 + 2X)m A - X(l + X)m 



(25) 



4(toa + m) 

The resulting numbers for Sa A vary between —6 ■ 10 -4 fm 3 and +4 ■ 10 -4 fm 3 for 
the ranges —0.8 < X, Y < 0.4, 4 < g± < 5 and 4.5 < g 2 < 9.5. In the absence of 
more stringent bounds on these parameters, we will set da A n {m A ) = and assign the 
theoretical predictions for the electric polarizabilies an error of ±2-10~ 4 fm 3 accordingly. 

Another contribution to the coefficients 5ai(X) and 8 fa (A) comes from loops involv- 
ing charged kaons [24]. Since we are working in SU(2), the kaons and etas are frozen 
out and effectively give some finite contact terms. To improve the estimate given in [1] , 
we include the average nucleon-A, E° mass splitting, 



A = ]- (m A + m E o - 2m) =216 MeV 
For the electric and magnetic polarizabilities, this leads to 



(26) 



5a* (m A ) = ^(D 2 + 3F 2 ) F 1 (M K , A), 



C 



5(3* (m A ) = ^(D 2 + 3F 2 ) F 2 (M Kl A), 



6"7T 

C 



5a«(m A ) = ^-(D-F) 2 F 1 (M K ,A) 



C 



6n 



8ft (m A ) = ^( D - F ) Fi{M K , A) 



Fi(Mk, A) = 
F 2 (M K ,A) = 



9A 



A 2 -Ml 



+ 



lOM^ - A 2 

-A 2 ) 3 / 2 
A 



(27) 



arccos 



M 



K 



VM 2 K - A' 



arccos 



M 



K 



3 - ) Notice that such an approach does not have a consistent chiral power counting as 
shown in [23] . It might be justified in the limit of an infinite number of colors where the 
nucleon and the A are degenerate in mass. 
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If one sets m\ = m^o = m (which is consistent to the order we are working) , one recovers 
eq.(14) of ref.[l]. With D = 0.8, F = 0.5 and M K = 495 MeV one finds 5a£(m A ) = 
1.31 • 10 -4 fm 3 and da^(rriA) = 0.13-10 -4 fm 3 . The corresponding numbers for the kaon 
contributions to the magnetic polarizabilities are a factor F2/F1 = 0.12 smaller. In the 
case of mass degeneracy for baryons this factor is exactly 1/10 [24]. The values based on 
(27) might, however, considerably overestimate the kaon loop contribution. Integrating 
e.g. the data from ref.[27] for — > K A, KTP , one gets a much smaller contribution 
since the typical cross sections are of the order of a few jubarn. This points towards the 
importance of a better understanding of SU(3) breaking effects. At present, we can not 
offer a solution to resolve this discrepancy. 4 ^ 

We now are in the position to give a prediction for the electromagnetic polariz- 
abilities. As our central (best) values we take the ones from [1]. The main sources of 
uncertainty stem from the choice of the value for g_&, the scale A and the kaon and A 
contributions to the various low-energy constants. Adding these in quadrature leads to 

a p = (10.5 ± 2.0) • 10" 4 fm 3 % = (3.5 ± 3.6) • 10" 4 fm 3 
a n = (13.4 ± 1.5) • 10" 4 fm 3 % = (7.8 ± 3.6) • 10" 4 fm 3 

which with the exception of j3 n agree with the empirical data [25] 

a p = (10.4 ± 0.6) • 10" 4 fm 3 , (3 P = (3.8 q= 0.6) • 10" 4 fm 3 
a n = (12.3 ±1.3) -10- 4 fm 3 , f3 n = (3.5 =f 1.3) • 10" 4 fm 3 . 

making use of the dispersion sum rules [12,26] (a + f3) p = (14.2 ± 0.3) ■ 10~ 4 fm 3 and 
(a + j3) n = (15.8 ± 0.5) • 10~ 4 fm 3 . 5 ^ Notice that we have added the systematic and 
statistical errors of the empirical determinations in quadrature. Clearly, an independent 
determination of the electric and magnetic nucleon polarizabilities would be needed to 
further tighten the empirical bounds on these fundamental quantities. This was also 
stressed in ref.[22]. It is worth to point out that the uncertainties given in (28) do not 
include effects of two (and higher) loops which start out at order q 5 . We do not expect 
these to alter the prediction for the electric polarizabilities significantly [1]. Such an 
investigation is underway but goes beyond the scope of this paper. Notice also that 
at present the theoretical uncertainties are larger than the experimental ones (if one 

4 ) We are grateful to Anatoly L'vov for drawing our attention to this problem. 

5 ) Notice that the uncertainty on the sum rule value for the neutron is presumably 
underestimated since one has to use deuteron data to extract the photon-neutron cross 
section. 
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imposes the sum rules for (a + ($)). That there is more spread in the empirical numbers 
when the dispersion sum rules are not used can e.g. be seen in the paper of Federspiel 
et al. in ref.[25]. 
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FIGURE CAPTIONS 

Fig.l Feynman rules for the 1/m suppressed propagators and vertices. Solid, dashed 
and wiggly lines refer to nucleons, pions and photons, in order. Pion momenta are 
denoted by (i = 1,2,3), nucleon momenta by £, p\, p%, and isospin indices by 

a, 6, a 

Fig. 2 The spin-averaged Compton amplitude for the proton normalized to one at ui = 0. 
(a) A(z) and A(() compared to the data for < oo < 150 MeV. (b) A(z) and A(Q 
compared to the quadratic approximation and to the data for 70 < u < 150 MeV. 

Fig. 3 The spin-averaged Compton amplitude for the neutron. 

Fig. 4 The real part of A p (uj) normalized to one at to = for uj > 140 MeV in comparison 
to the data. 

Fig. 5 Dependence of the electromagnetic polarizabilities on ci for 5ai(m/\) = 5f3i(m&) = 
(i = p, n) in units of 10~ 4 fm 3 . 

Fig. 6 Dependence of the electromagnetic polarizabilities on A for 5ai(\) = 5j3i(\) = 
(i = p,n) in units of 10 fm 3 . The upper/lower rim of the corresponding bands 
refers to gA = 1.328 / 1.26, in order. 

Fig.7 5/3 A (#i, Y) in units of 10" 4 fm 3 for 3.8 < g x < 5.3 and -0.8 < Y < 0.3. 
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